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ABSTRACT 

This  study  evaluates  the  intraseasonal  variahility  associated  with  summer  precipitation  over  South  America 
in  14  coupled  general  circulation  models  (GCMs)  participating  in  the  Intergovernmental  Panel  on  Climate 
Change  (IPCC)  Fourth  Assessment  Report  (AR4).  Eight  years  of  each  modePs  twentieth-century  climate 
simulation  are  analyzed.  Two  dominant  intraseasonal  hands  associated  with  summer  precipitation  over  South 
America  are  focused  on:  the  40-  and  the  22-day  band.  The  resulis  show  that  in  the  souihern  summer 
(November-April),  most  of  the  models  underestimate  seasonal  mean  precipitation  over  central-east  Brazil, 
northeast  Brazil,  and  the  South  Atlantic  convergence  zone  (SACZ),  while  the  Atlantic  intertropical  con¬ 
vergence  zone  (ITCZ)  is  shifted  southward  of  its  observed  position.  Most  of  the  models  capture  both  the 
40-  and  22-day  band  around  Uruguay,  but  with  less  frequent  active  episodes  than  observed.  The  models  also 
tend  to  underestimate  the  total  intraseasonal  (10-90  day),  the  40-,  and  the  22-day  band  variances.  For  the 
40-day  band,  10  of  the  14  models  simulate  to  some  extent  the  3-cell  pattern  around  South  America,  and 
6  models  reproduce  its  teleconnection  with  precipitation  in  the  south-central  Pacific,  but  only  I  model  sim¬ 
ulates  the  teleconnection  with  the  MJO  in  the  equatorial  Pacific,  and  only  3  models  capture  its  northward 
propagation  from  50°  to  32°S.  For  the  7  models  with  three-dimensional  data  available,  only  1  model  repro¬ 
duces  well  the  deep  baroclinic  vertical  structure  of  the  40-day  band.  For  the  22-day  band,  only  6  of  the 
14  models  capture  its  northward  propagation  from  the  SACZ  to  the  Atlantic  ITCZ.  It  is  found  that  models 
with  some  form  of  moisture  convective  trigger  tend  to  produce  large  variances  for  the  intraseasonal  bands. 


1.  Introduction 

The  climate  of  tropieal  South  America  is  characterized 
by  a  pronounecd  summer  monsoon,  whieh  is  often  re¬ 
ferred  to  as  the  South  Ameriean  monsoon  system  (SAMS; 
Kousky  1988;  Horel  et  al.  1989;  Lenters  and  Cook  1995; 
Zhou  and  Lau  1998;  see  reviews  by  Nogues-Paegle  et  al. 
2002;  Vera  et  al.  2006a).  The  summer  precipitation  over 
South  America  has  strong  intraseasonal  variability  with 
the  leading  pattern  of  deep  convection  showing  a  seesaw 
between  the  South  Atlantie  convergence  zone  (SACZ) 
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and  the  subtropical  plains  of  South  America  (Nogues- 
Paegle  and  Mo  1997)  Further  studies  show  that  this 
seesaw  pattern  is  part  of  a  mueh  larger  Rossby  wave 
train  structure  that  include  alternating  eenters  of  negative 
and  positive  streamfunction,  geopotential  height,  and 
temperature  anomalies  in  the  southern  portion  of  the 
eontinent,  and  farther  upstream  in  the  southern  Paeifie 
(Liebmann  et  al.  1999, 2004;  Paegle  et  al.  2000;  Jones  and 
Carvalho  2002;  Diaz  and  Aceituno  2003;  Carvalho  et  al. 
2004).  Using  singular  speetrum  analysis,  Paegle  et  al. 
(2(XX))  found  that  this  seesaw  pattern  is  dominated  by  two 
frequency  bands:  a  band  with  a  period  of  about  36-40 
days  (hereafter  the  40-day  band)  and  a  band  with  a  period 
of  about  22-28  days  (hereafter  the  22-day  band).  Both 
bands  are  linked  to  tropical  convection.  The  40-day  band 
is  related  to  the  Madden-Julian  oscillation  (MJO)  in  the 
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tropics  while  the  22-day  band  is  connected  to  a  tropical 
mode  at  the  corresponding  frequency  band.  When  the 
SACZ  is  enhanced,  these  two  bands  become  meridionally 
aligned  locally  and  such  episodes  are  characterized  by  a 
wave  train  propagating  northward  from  southern  South 
America  toward  the  tropics.  These  intrascasonal  bands 
are  responsible  for  alternating  wet  and  dry  episodes  over 
the  SAMS  region.  Few  studies,  however,  have  evaluated 
the  simulations  of  the  intrascasonal  variability  of  the 
summer  precipitation  over  South  America  by  general 
circulation  models  (GCMs).  Misra  (2005)  examined  the 
simulation  by  one  atmospheric  GCM  and  found  the  in- 
traseasonal  variability  to  be  inadequately  represented. 
Furthermore,  downsealing  the  GCM  results  to  a  regional 
model  did  not  improve  the  variability. 

Reeently,  in  preparation  for  the  Intergovernmental 
Panel  on  Climate  Change  (IPCC)  Fourth  Assessment 
Report  (AR4),  more  than  a  dozen  international  climate 
modeling  centers  conducted  a  comprehensive  set  of  long¬ 
term  simulations  for  both  the  twentieth-century’s  climate 
and  different  climate  change  scenarios  in  the  twenty-first 
century,  which  constitutes  the  third  phase  of  the  World 
Climate  Research  Programme  (WCRP)  Coupled  Model 
Intercomparison  Project  (CMIP3;  Mechl  et  al.  2007). 
This  is  an  unprecedented,  comprehensive  coordinated  set 
of  global  coupled  climate  experiments  for  the  twentieth 
and  twnty-first  century.  Before  conducting  the  extended 
simulations,  many  of  the  modeling  centers  applied  an 
overhaul  to  their  physical  schemes  to  incorporate  state- 
of-the-art  research  results.  For  example,  almost  all 
modeling  centers  have  implemented  prognostic  cloud 
microphysics  schemes  in  their  models,  some  have  added 
a  moisture  trigger  to  their  deep  convection  schemes,  and 
some  now  take  into  account  convective  momentum 
transport.  Moreover,  many  modeling  centers  increased 
their  models’  horizontal  and  vertical  resolutions  and 
some  conducted  experiments  with  different  resolutions. 

The  purpose  of  this  study  is  to  evaluate  the  intra- 
seasonal  variability  of  precipitation  associated  with  the 
summer  precipitation  over  South  America  in  14  IPCC 
AR4  coupled  GCMs,  with  emphasis  on  the  40-  and  the 
22-day  bands.  While  there  has  been  some  analysis  of 
seasonal  means  of  the  IPCC  runs  in  this  region  (e.g., 
Vera  ct  al.  2006b),  intraseasonal  variability  has  not  been 
studied.  The  models  and  validation  datasets  used  in  this 
study  are  described  in  section  2.  The  diagnostic  methods 
and  results  are  described  in  section  3.  A  summary  and 
discussion  are  given  in  section  4. 

2.  Models  and  validation  datasets 

This  analysis  is  based  on  8  years  of  the  climate  of  the 
twentieth  century  (20C3M)  simulations  from  14  coupled 


GCMs.  Table  1  shows  the  model  names  and  acronyms, 
their  horizontal  and  vertical  resolutions,  and  brief  de¬ 
scriptions  of  their  deep  convection  schemes.  For  each 
model  we  use  8  yr  of  daily  mean  surface  precipitation. 
Three-dimensional  data  arc  available  for  7  of  the  14 
models,  for  which  we  analyzed  upper-air  winds,  tem¬ 
perature,  and  specific  humidity. 

The  model  simulations  are  validated  using  the  Global 
Precipitation  Climatology  Project  (GPCP)  version  2 
precipitation  (Huffman  et  al.  2001).  We  use  8  yr  (1997- 
2004)  of  daily  data  with  a  horizontal  resolution  of 
1°  X  1°.  We  also  use  8  yr  (1997-2004)  of  daily  National 
Centers  for  Environmental  Prediction-National  Center 
for  Atmospheric  Research  (NCEP-NCAR)  reanalysis 
data  (Reanalysis  I;  Kalnay  et  al.  1996),  for  which  we 
analyzed  upper-air  winds,  temperature,  and  specific 
humidity. 

Total  intrascasonal  (periods  of  10-90  days)  anomalies 
were  obtained  by  applying  a  365-point  10-90-day  Lanezos 
filter  (Duchan  1979).  Because  the  Lanezos  filter  is  non- 
recursivc,  182  days  of  data  were  lost  at  each  end  of  the 
time  series  (364  days  in  total).  The  dominant  intra¬ 
seasonal  bands  are  determined  using  wavelet  spectrum 
because  they  are  active  mainly  during  the  southern 
summer.  Wavelet  spectrum  is  a  powerful  tool  for  an¬ 
alyzing  multiscale,  nonstationary  processes,  and  can 
simultaneously  determine  both  the  dominant  bands  of 
variability  and  how  those  bands  vary  in  time  (e.g.,  Mak 
1995;  Torrence  and  Compo  1998).  We  utilize  the  wavelet 
analysis  program  developed  by  Torrence  and  Compo 
(1998)  and  use  the  Morlet  wavelet  as  the  mother  wavelet. 
The  40-day  band  is  defined  as  precipitation  variability  in 
the  period  range  of  30-60  days,  and  was  obtained  by 
applying  a  365-point  30-60-day  Lanezos  filter.  Similarly, 
the  22-day  band  is  defined  as  precipitation  variability  in 
the  period  range  of  20-30  days,  again  using  a  365-point 
Lanezos  filter.  We  also  tested  the  Murakami  (1979)  filter 
with  similar  results. 

3.  Results 

a.  Southern  summer  (November- April)  seasonal 
mean  precipitation 

Previous  observational  studies  indicate  that  the  intra¬ 
seasonal  variance  of  precipitation  is  highly  correlated 
with  time-mean  precipitation  (e.g.,  Wheeler  and  Kiladis 
1999).  That  is,  areas  with  abundant  mean  precipitation 
tend  to  be  characterized  by  large  intraseasonal  variabil¬ 
ity.  Therefore,  we  first  look  at  the  horizontal  distribution 
of  southern  summer  (November- April)  seasonal  mean 
precipitation  [Fig.  1;  see  also  Vera  et  al.  (2006b)  for  an 
evaluation  of  3-month  season  climatologies  of  the  IPCC 
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Table  1.  List  of  models  that  participate  in  this  study. 


Modeling  groups 

IPCCID  (label  in 
figures) 

Grid 

type/resolution/model  top 

Deep  convection 
sehcme/modification 

Downdrafts* 

SaUC/Meso 

Closure/ 

trigger 

NOAA/GFDL 

GFDL-CM2.0 

(GFDL2.0) 

Grid  point 

144/90/L24  3  mb 

Moorthi  and  Suarez 
(1992)/Tokioka 
et  al.  (1988) 

N/N/N 

CAPE/threshold 

NOAA/GFDL 

GFDL-CM2.1 

(GFDL2.1) 

Grid  point 

I44/90/L24  3  mb 

Moorthi  and  Suarez 
(1992)/Tokioka 
et  al.  (1988) 

N/N/N 

CAPE/threshold 

NCAR 

CCSM3  (CCSM3) 

Spectral  T85/L26  2.2  mb 

Zhang  and  MeFarlane 
(1995) 

Y/N/N 

CAPE 

NCAR 

PCM  (PCM) 

Spectral  T42/L26  2.2  mb 

Zhang  and  MeFarlane 
(1995) 

Y/N/N 

CAPE 

NASA/GISS 

GISS-AOM 

(GISS-AOM) 

Grid  point  90/60/Ll 2 

Russell  et  al.  (1995) 

N/N/N 

CAPE 

NASA/GISS 

GlSS-ER  (GISS-ER) 

Grid  point 

72/46/L20  0.1  mb 

Del  Genio  and  Yao 
(1993) 

Y/N/N 

Cloud  base 
buoyancy 

Center  for  Climate 
System  Research, 
National  Institute 

MlROC3.2(hires) 

(MIROC-hircs) 

Spectral  T106/L56 

Pan  and  Randall 
(1998)/Emori 
et  al.  (2001) 

Y/N/N 

CAPE/relative 

humidity 

for  Environmental 

Studies,  and  Frontier 

Research  Center 
for  Global  Change 

Same  as  above 

M I  ROC3.2(  medres) 
(MIROC-medres) 

Spectral  T42/L20  30  km 

Pan  and  Randall 
(1998)/Emori 
et  al.  (2001) 

Y/N/N 

CAPE/relative 

humidity 

MRI 

MRI-CGCM2.3.2 

(MRI) 

Spectral  T42/L30  0.4  mb 

Pan  and  Randall 
(1998) 

Y/N/N 

CAPE 

Canadian  Centre  for 
Climate  Modeling 

CGCM3.1-T47 

(CGCM) 

Spectral  T47*L32  1  mb 

Zhang  &  MeFarlane 
(1995) 

Y/N/N 

CAPE 

and  Analysis 

MPI 

ECHAM5/MPI-OM 

(MPI) 

Spectral  T63/L31  10  mb 

Tiedtke  (1989)/ 
Nordeng  (1994) 

Y/N/N 

CAPE/moisture 

eonvergenee 

IPSL 

IPSL-CM4(IPSL) 

Grid  point  96/72/L19 

Emanuel  (1991) 

Y/Y/N 

CAPE 

Meteo-France/CN  R  M 

CNRM-CM3 

(CNRM) 

Spectral  T63/L45 

0.05  mb 

Bougeault  (1985) 

N/N/N 

Kuo 

CSIRO  Atmospheric 
Research 

CSIRO  Mk3.0 
(CSIRO) 

Spectral  T63/L18  4  mb 

Gregory  and 

Rowntree  (1990) 

Y/N/N 

Cloud  base 
buoyancy 

*  For  downdrafts,  SC  means  saturated  convective  downdrafts,  UC  means  unsaturated  convective  downdrafts,  and  Meso  means  mesoscalc 
downdrafts;  indicated  by  yes  (Y)  or  no  (N). 


runs].  The  observed  large-scale  preeipilalion  (Fig.  la) 
pattern  is  one  of  intense  precipitation  over  the  Amazon 
basin,  an  eastern  Paeifie  intertropieal  eonvergenee  zone 
(ITCZ),  an  Atlantie  ITCZ,  plus  a  band  of  enhaneed 
precipitation  that  extends  to  the  southeast  from  the 
maximum  in  the  Amazon,  known  as  the  SACZ  (e.g., 
Kodama  1992). 

Most  of  the  models  underestimate  preeipitation  over 
the  Amazon  basin.  Only  a  few  models  produee  magni¬ 
tude  elose  to  that  observed  (i.e.,  the  Model  for  Inter- 
diseiplinary  Researeh  on  Climate  3.2,  high-resolution 
version  [MIROC3.2(hires)]  and  the  Parallel  Climate 
Model  (PCM)|.  The  maximum  is  shifted  to  the  east  in 
three  models  [i.e.,  the  Geophysieal  Fluid  Dynamics 
Laboratory  model  version  2.0  (GFDL2.0),  the  Goddard 


Institute  for  Space  Studies  Atmosphere-Oeean  Model 
(GISS-AOM),  and  the  Goddard  Institute  for  Spaee 
Studies  Model  E-R  (GISS-ER)].  There  are  loeal  max¬ 
ima  over  the  Andes  Mountains  in  eight  models  (i.e., 
PCM,  GISS-AOM,  the  Model  for  Interdiseiplinary 
Researeh  on  Climate  3.2,  medium-resolution  version 
[MIROC3.2(medres)],  MIROC3.2(hires),  the  Meteo- 
rologieal  Researeh  Institute  Coupled  General  Cireula- 
tion  Model,  version  2.3.2a  (MRI  CGCM),  L’Institut 
Pierre-Simon  Laplaee  (IPSL),  Centre  National  de  Re- 
eherehes  Meteorologiques  (CNRM))  that  do  not  exist 
in  GPCP.  However,  it  is  important  to  note  that  some 
other  preeipitation  analyses  for  South  Ameriea  depiet  a 
preeipitation  maximum  along  the  tropieal  Andes  (e.g., 
Hoffman  1975).  The  eastern  Paeifie  ITCZ  is  shifted  south 
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Fig.  1  Southern  summer  (November-April)  seasonal  mean  precipitation  for  observation  and  14  IPCC  AR4  models. 
The  first  contour  is  4  mm  day  ’  and  the  contour  interval  is  2  mm  day 
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of  the  equator  in  five  models  [i.e.,  the  Community  Cli¬ 
mate  System  Model,  version  3  (CCSM3),  PCM,  GISS- 
ER,  IPSL,  Commonwealth  Scientific  and  Industrial 
Research  Organisation  (CSIRO)]  and  there  is  a  double- 
ITCZ  pattern  in  the  eastern  Pacific  in  six  models  [i.e., 
GFDL2.0,  GFDL2.1,  GISS-AOM,  MIROC3.2(medres), 
MIROC3.2(hires),  CNRM;  see  also  Lin  2007].  The  At¬ 
lantic  ITCZ  is  too  far  south  in  almost  all  models,  and  two 
models  (i.e.,  GFDL2.0  and  GFDL2.1)  show  a  double- 
ITCZ  pattern  in  the  tropical  Atlantic.  Finally,  simulated 
precipitation  in  the  S ACZ  is  almost  always  too  weak,  and 
in  the  seven  models  that  do  contain  an  S  ACZ  signature  it 
is  shifted  northward  with  respect  to  observations  [i.e., 
GFDL2.0,  GFDL2.1,  PCM,  MIROC3.2(hires),  MRI,  the 
Max  Planck  Institute  (MPI),  and  CNRMj. 

As  will  be  shown  shortly,  the  largest  intraseasonal 
variability  associated  with  summer  precipitation  over 
South  America  is  concentrated  in  a  meridional  belt  be¬ 
tween  30°  and  60°W  (roughly  the  eastern  continent  and 
the  western  Atlantic  Ocean).  Therefore,  we  conduct  a 
more  quantitative  evaluation  of  the  seasonal  mean 
precipitation  averaged  over  these  longitudes  (Fig.  2). 
Observations  reveal  two  local  maxima:  one  at  2°S  cor¬ 
responding  to  the  Amazon  precipitation  and  Atlantic 
ITCZ,  and  a  secondary  peak  at  30°S  corresponding  to 
the  SACZ.  Almost  all  of  the  models  show  only  one 
maximum.  In  total  11  models  have  their  maximum 
shifted  southward  compared  to  observed,  to  10°S  [i.e., 
GFDL2.0,  CCSM3,  GISS-AOM,  MIROC3.2(hires), 
MRI,  CGCM,  MPI,  IPSL,  CSIRO]  or  15°S  (i.e.,  PCM 
and  CNRM),  which  is  associated  with  overly  weak 
Amazon  precipitation,  and/or  the  southward  shift  of 
Amazon  precipitation/ Atlantic  ITCZ  in  those  models. 
All  models  underestimate  the  precipitation  at  30°S,  which 
is  often  associated  with  a  too  weak  SACZ  extension  into 
the  Atlantic.  For  the  region  between  10°  and  25°S,  9  of  the 
14  models  produce  quite  reasonable  precipitation  [i.e., 
GFDL2.0,  CCSM3,  GISS-AOM,  MIROC3.2(medres), 
MIROC3.2(hires),  MRI,  CGCM,  MPI,  IPSL,  CSIRO], 
while  two  models  overestimate  precipitation  (i.e.,  PCM 
and  CNRM)  and  three  models  underestimate  it  (i.e., 
GFDL2.1,  GISS-ER,  and  IPSL). 

b.  Total  intraseasonal  (10-90  day)  variance 

Figure  3  shows  the  horizontal  distribution  of  the 
standard  deviation  of  total  intraseasonal  (10-90  day) 
precipitation  anomaly  during  the  November- April  sea¬ 
son.  In  observations  (Fig.  3a),  the  intraseasonal  variance 
does  not  completely  follow  that  of  seasonal  mean  pre¬ 
cipitation  (Fig.  la),  but  is  concentrated  from  approxi¬ 
mately  10°N-40°S  between  30°  and  60°W.  There  are  three 
local  maxima:  over  the  Amazon  River  mouth,  over  the 
Atlantic  extension  of  the  SACZ,  and  over  southeast 
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Latitude  (deg) 


Fig.  2.  Meridional  profile  of  soulhern  summer  (November- 
April)  seasonal  mean  precipitation  averaged  between  30®  and  60®W 
for  the  observations  and  the  14  models. 

Brazil/Uruguay.  These  are  consistent  with  the  results 
of  Liebmann  et  al.  (1999).  The  mismatch  between  the 
seasonal  mean  precipitation  and  total  intraseasonal 
variance  suggests  that  the  intraseasonal  variability  is 
more  than  simply  noise  around  the  seasonal  mean,  but  is 
caused  by  mechanisms  that  vary  from  those  related  to 
seasonal  mean  precipitation.  Therefore,  it  is  of  interest 
as  to  whether  if  the  models  are  able  to  reproduce  this 
mismatch.  The  model  variances  show  two  eharaeteris- 
tics.  First,  in  8  of  the  14  models  the  distribution  of  in¬ 
traseasonal  variance  does  not  completely  follow  that  of 
the  seasonal  mean  precipitation  [i.e.,  GFDL2.0,  GFDL2.1, 
MIROC3.2(hires),  MRI,  CGCM,  MPI,  IPSL,  and  CSIRO]. 
In  three  models  the  intraseasonal  variance  follows  the 
mean  precipitation  [i.e.,  PCM,  MIROC3.2(medres),  and 
CNRM],  and  in  three  models  the  intraseasonal  variance 
is  too  small  (i.e.,  CCSM3,  GISS-AOM,  and  GISS-ER). 
Second,  the  models  tend  to  produce  their  maximum 
variance  over  their  SACZ,  but  fail  to  produce  the  maxima 
over  the  Amazon  River  mouth  or  Uruguay. 

To  provide  a  more  quantitative  evaluation  of  the 
model  simulations.  Fig.  4  shows  the  meridional  profile  of 
total  intraseasonal  (10-90  day)  variance  of  precipitation 
averaged  between  30°  and  60°W.  The  observed  variance 
shows  three  peaks  at  2°,  17°,  and  32°S.  All  models  un¬ 
derestimate  the  variance  around  2°  and  32°S.  Only  a  few 
of  the  models  produce  any  sort  of  a  peak  at  all  in  those 
regions.  For  the  region  between  10°  and  25°S,  six  models 
simulate  nearly  realistic  or  overly  large  variance  [i.e., 
MIROC3.2(medres),  MIROC3.2(hires),  MPI,  CNRM, 
CSIRO,  and  GFDL2.0].  The  other  eight  models  un¬ 
derestimate  variance,  although  six  of  the  eight  display 
reasonable  seasonal  mean  precipitation  in  this  region 
(Fig.  2).  Interestingly,  the  six  models  simulating  nearly 
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Fig.  3.  Horizontal  distribution  of  the  standard  deviation  of  total  intraseasonal  (10-90  day)  precipitation  anomaly 
during  northern  summer  (November-April).  The  first  contour  is  3  mm  day"  ‘  and  the  contour  interval  is  1  mmday  *. 
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Fig.  4.  Meridional  profile  of  the  total  intraseasonal  (10-90  day) 
variance  of  precipitation  anomaly  averaged  between  30°  and  60°W. 

realistic  or  overly  large  variance  are  the  same  models 
that  contain  large  variances  for  the  convectively  coupled 
equatorial  waves  (Lin  et  al.  2006).  A  common  charac¬ 
teristic  of  these  models  is  that  there  is  some  form  of 
moisture  trigger  of  their  convection  scheme,  suggesting 
that  a  moisture  trigger  for  deep  convection  may  improve 
the  simulation  of  intraseasonal  variability  associated 
with  summer  precipitation  over  South  America. 

c.  The  dominant  intraseasonal  bands 

Figure  5  shows  the  wavelet  spectrum  of  precipita¬ 
tion  averaged  between  30°-35°S,  50°“60°W  (around 
Uruguay)  for  observations  and  the  14  IPCC  models. 
TTie  observed  spectrum  (Fig.  5a)  demonstrates  two  dom¬ 
inant  intraseasonal  bands,  a  30-60-day  band  (the  so-called 
40-day  band)  and  a  20-30-day  band  (the  so-callcd  22-day 
band).  Most  of  the  models  capture  both  bands,  although 
the  model  variances  are  generally  smaller  than  the 
observed  variances.  The  models  also  tend  to  produce 
fewer  active  episodes.  Only  two  models  (i.e.,  CGCM 
and  CSIRO)  produce  frequent  active  episodes  in  both 
bands.  It  is  important  to  note  that  many  models  have 
excessively  large  power  between  60  and  100  days.  This 
suggests  that  the  summer  precipitation  in  the  models 
has  larger  persistence  than  is  observed.  Lin  et  al.  (2(X)6) 
found  a  similar  problem  associated  with  the  tropical 
oceanic  precipitation  in  the  models,  and  hypothesized 
that  it  is  caused  by  the  erroneous  representation  of  self¬ 
suppression  processes  in  deep  convection  in  the  modeLs 
moisture  physics. 

d.  The  40-day  hand 

Next  we  focus  on  the  40-day  band.  Figure  6  shows  the 
meridional  profile  of  the  40-day  band  variance  averaged 


between  30°  and  60°W.  The  observed  variance  in  the 
40-day  band  is  similar  to  the  total  intraseasonal  variance 
in  that  there  are  peaks  at  2°,  17°,  and  32°S.  At  40  days, 
however,  the  maximum  at  17°S  is  larger  than  that  at 
32°S,  while  for  the  total  intraseasonal  variance  (Fig.  4) 
the  17°S  peak  is  relatively  small,  and  is  about  the  same  as 
that  at  32°S,  All  models  underestimate  the  variance  near 
2°  and  32°S.  For  the  region  between  10°  and  25°S,  the 
six  models  producing  realistic  or  excessive  total  intra- 
seasonal  variance  produce  the  40-day  band  variance  that 
is  between  the  observed  value  and  half  the  observed 
value,  and  the  wavelet  analysis  (Fig,  5)  suggests  that 
their  intraseasonal  variance  is  concentrated  more  in  the 
lower-frequency  band.  The  other  eight  models  produce 
the  40-day  variance  that  is  less  than  half  of  the  observed 
value,  although  six  of  them  display  reasonable  seasonal 
mean  precipitation  in  this  region  (Fig.  2).  Possible  rea¬ 
sons  of  this  will  be  discussed  in  section  4.  It  is  important 
to  note  that  although  the  models  do  not  simulate  the 
right  intensity  of  variance,  some  are  able  to  simulate  the 
position  of  its  peaks  (e.g.,  GFDL2.0  and  MPI). 

Figure  7  shows  the  lag  correlation  of  40-day  bandpass- 
filtered  precipitation  at  30°S,  55°W  with  the  40-day 
precipitation  averaged  from  50°  to  60°W.  Shading  de¬ 
notes  the  regions  above  the  95%  confidence  level.  The 
observations  (Fig.  7a),  as  expected,  show  a  40-day  pe¬ 
riod  oscillation  at  the  latitude  of  the  base  grid  point. 
Zonal  anomalies  of  the  opposite  sign  centered  at  about 
15°S  slightly  lead  those  at  the  base  grid  point,  resulting  in 
a  dipole  pattern.  Although  the  choice  of  base  grid  point 
makes  the  anomalies  appear  strongest  at  30°  and  15°S, 
the  figure  shows  that  the  dipole  actually  propagates 
northward,  with  the  opposite-signed  anomalies  that 
are  evident  at  the  latitude  of  the  base  grid  point,  but 
20  days  prior,  moving  northward  to  become  the  near- 
simultaneous  antinodc  some  20°  to  the  north.  The  ob¬ 
served  northward  propagation  of  the  dipole  is  consistent 
with  the  results  of  Nogues-Paegle  and  Mo  (1997)  and 
Diaz  and  Accituno  (2003). 

Only  2  of  the  14  models  capture  both  the  north¬ 
ward  propagation  and  the  dipole  [i.e.,  GISS-AOM 
and  MIROC3.2(hircs)].  One  model  captures  only  the 
northward  propagation  (MPI),  while  six  models  cap¬ 
ture  only  the  dipole  [i.e.,  GFDL2.1,  CCSM3,  PCM, 
MIROC3.2(medres),  MRI,  and  CGCM].  The  other  five 
models  lack  either  of  these  features.  It  is  important  to 
note  that  the  reason  some  models  (e.g.,  GFDL2.0)  do  not 
show  the  dipole  structure  linking  the  subtropics  to  the 
tropics  is  because  the  tropical  center  in  the  model  simu¬ 
lations  is  outside  of  the  band  50°-60°W.  For  example,  it 
will  be  shown  (Fig,  11b)  that  such  a  dipolar  structure  is 
present  in  GFDL2.0,  even  though  it  is  not  evident  in 
Fig.  7b. 
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Fig.  5.  Wavelet  spectrum  of  precipitation  averaged  between  30°-35°S»  50°-60°W.  Only  spectral  peaks  above  the  95% 

confidence  level  are  plotted. 
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Fig.  6.  As  in  Fig.  4,  but  for  the  variance  of  the  40-day  band. 

Next  wc  examine  the  vertieal  struetures  of  the  models 
in  the  40-day  band.  Figure  8  shows  the  lag  correlation  of 
temperature  averaged  between  20'^-30°S  and  50'^-60°W 
versus  the  40-day  band  preeipitation  anomaly  at  the 
same  loeation  for  observation  and  the  seven  models  with 
three-dimensional  data  available.  Note  that  for  four 
models  the  3D  data  extends  to  only  200  mb.  In  obser¬ 
vations,  the  40-day  band  displays  a  deep  warm  core 
between  the  surfaee  and  200  mb  and  a  eold  eore  above 
250  mb  during  the  eonveetive  phase.  Six  of  the  seven 
models  show  a  signifieant  warm  temperature  anomaly, 
but  often  with  a  large  southward  phase  tilt  with  height. 

There  is  also  a  signifieant  bias  in  the  geopotential 
height  structure  in  many  models  (Fig.  9).  The  observed 
geopotential  height  displays  a  deep  baroelinic  structure, 
with  a  positive  anomaly  extending  from  the  tropopause 
to  750  mb  and  a  negative  anomaly  from  750  mb  to  the 
surface  during  the  convective  phase  (Fig.  9a).  Only  one 
model  (i.e.,  GFDL2.0)  reproduces  the  deep  baroelinic 
structure.  In  the  other  six  models  the  negative  anomaly 
extends  too  high  into  the  middlc/upper  troposphere, 
indicating  a  more  barotropic  structure. 

Figure  10  shows  the  vertical  structure  of  divergence. 
The  observed  divergence  displays  a  two-layer  structure 
during  the  preeipitating  phase,  with  convergence  from 
the  surfaee  to  450  mb,  and  divergenee  above  450  mb 
(Fig.  lOa).  All  but  one  model  (i.e.,  MPI)  reproduce 
fairly  well  the  two-layer  strueture,  although  in  GFDL2.1 
(Fig.  lOe)  the  eonvergence  layer  is  too  deep,  extending 
from  the  surface  to  350  mb.  Previous  studies  (e.g.,  Paegle 
ct  al.  2000)  show  that  precipitation  variability  in  the 
30-60-day  band  observed  at  the  region  around  30'^S, 
55°W  is  assoeiated  with  the  aetivity  of  Rossby  wave 
trains  propagating  into  the  region  from  the  South  Pacifie. 
Therefore,  it  seems  that  biases  assoeiated  with  tempera¬ 
ture,  geopotential  height,  and  divergence  are  related  to 


modeling  deficiencies  in  reproducing  the  features  asso¬ 
ciated  with  the  Rossby  wave  trains. 

Next  we  look  at  the  teleconnection  pattern  associated 
with  the  40-day  band.  Figure  11  shows  the  linear  eorre- 
lation  of  the  40-day  band  precipitation  anomaly  versus 
itself  averaged  between  25®-35°S  and  30°-60°W.  In  ob¬ 
servations  (Fig.  11a),  there  is  a  three-cell  pattern  around 
South  America  with  a  positive  precipitation  anomaly 
over  Uruguay  and  negative  anomalies  over  the  SACZ 
and  the  South  Pacific  around  50°S,  280°E,  which  are  all 
statistically  significant  above  the  95%  confidence  level. 
This  three-cell  pattern  has  been  found  in  previous  ob¬ 
servational  studies  using  OLR  (Carvalho  et  al.  2004, 
see  their  Fig.  8c)  and  upper-air  geopotential  height, 
streamfunction,  and  winds  (Liebmann  et  al.  1999,  2004; 
Diaz  and  Aceituno  2003;  Carvalho  et  al.  2004).  At  the 
same  time,  there  is  a  dipole  over  the  tropical  Pacific  with 
a  negative  anomaly  over  the  eentral  Paeifie  and  positive 
anomaly  over  the  Maritime  Continent/western  Paeifie. 
These  are  consistent  with  the  results  of  Paegle  et  al. 
(2000,  see  their  Fig.  6d),  and  they  demonstrated  that  the 
dipole  over  tropieal  Paeifie  is  assoeiated  with  the  MJO. 
There  is  also  a  positive  anomaly  over  south  central  Pa¬ 
cific  around  20*^5, 200°E  with  a  SPCZ  developed  farther 
east  of  its  climatological  position,  which  is  consistent 
with  previous  work  (e.g.,  Nogues-Paegle  and  Mo  1997). 
In  total,  10  of  the  14  models  simulate  to  some  extent 
the  three-cell  pattern  around  the  South  America  [i.e., 
GFDL2.0,  GFDL2.1,  GISS-AOM,  MIROC3.2(mcdrcs), 
MRI,  CGCM,  MPI,  IPSL,  CNRM,  and  CSIRO]. 
However,  only  one  model  (i.e.,  GFDL2.0)  simulates  the 
MJO  dipole  over  the  tropical  Pacific.  Six  models  [i.e., 
GFDL2.1,  PCM,  GISS-ER,  MIROC3.2(hircs),  CGCM, 
MPI,  and  IPSL]  produce  statistically  a  significant  posi¬ 
tive  anomaly  in  the  south-central  Pacific  around  20°S, 
200"E. 

To  summarize,  all  models  substantially  underestimate 
the  40-day  band  variance  over  north  Brazil  and  Uru¬ 
guay,  while  about  half  of  the  models  simulate  nearly 
realistic  variance  over  the  SACZ.  In  total,  10  of  the 
14  models  simulate  to  some  extent  the  three-cell  pattern 
around  the  South  America,  with  6  models  reproducing 
its  telcconnection  with  precipitation  in  the  south-central 
Pacific.  However,  only  one  model  simulates  the  tele¬ 
connection  with  the  MJO  in  equatorial  Pacific,  and  only 
three  models  capture  its  northward  propagation  from 
50°  to  32°S.  Of  the  seven  models  with  three-dimensional 
data  available,  only  one  reproduces  well  the  deep 
baroelinic  vertical  structure  of  the  40-day  band. 

e.  The  22-day  hand 

Figure  12  shows  the  meridional  profile  of  the  22-day 
band  precipitation  variance  averaged  between  30°  and 
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7.  Lag  correlation  of  the  40-day  band  precipitation  anomaly  averaged  between  50°  and  60°W  with  respect  to  itself 
at  30°S,  55°W.  Shading  denotes  the  regions  where  the  lag  correlation  is  above  the  95%  confidence  level. 
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Fig.  7.  (Continued) 
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Fig.  8.  Lag  correlation  of  temperature  averaged  between  20®-30°S  and  5n*’-60°W  vs  the  40-day  band  preeipitation 
anomaly  at  the  same  location  for  observation  (NCEP-NCAR  reanalysis)  and  the  seven  models.  Shading  denotes 
the  area  where  eorrelation  is  above  the  95%  eonfidenee  level,  with  dark  (light)  shading  for  positive  (negative) 
eorrelation. 
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Fig.  9.  As  in  Fig.  8,  but  for  geof>otential  height. 
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Fig.  10.  As  in  Fig.  8,  but  for  divergence. 
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Fig.  11.  Linear  correlation  of  the  40-day  band  precipitation  anomaly  vs  itself  averaged  between  25®-35°S  and 
50°-60®W.  Shading  denotes  the  area  where  correlation  is  above  the  95%  confidence  level,  with  dark  (light)  shading 
for  positive  (negative)  correlation. 
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Fig.  12.  As  in  Fig.  4,  but  for  the  variance  of  the  22-day  band. 

60°W.  The  observed  profile  of  the  22-day  band  vari- 
anee  is  different  from  those  of  the  total  intraseasonal 
(10-90  day)  varianee  (Fig.  4)  and  the  40-day  band  vari- 
anee  (Fig.  6),  both  of  whieh  display  three  loeal  maxima 
with  the  primary  maximum  at  2°S.  The  22-day  band,  on 
the  other  hand,  shows  only  two  maxima  at  2°S  and  32°S 
with  the  later  having  slightly  larger  magnitude.  In  total, 
12  of  the  14  models  underestimate  the  varianee  around 
2°S,  and  all  models  underestimate  the  varianee  around 
32°S.  Between  10°  and  25°S,  six  models  simulate  realistie 
or  overly  large  varianee  [i.e.,  MIROC3.2(hires),  MPI, 
CNRM,  GFDL2.0,  MIROC3.2(medres),  and  CSIRO]. 
Again,  these  are  those  models  produeing  large  varianees 
for  the  eonveetively  eoupled  equatorial  waves  (Lin  et  al. 
2006). 

Figure  13  shows  the  lag-eorrelation  of  the  22-day  band 
preeipitation  anomaly  averaged  between  30°  and  60°W 
with  respeet  to  the  22-day  band  preeipitation  anomaly  at 
30°S,  55°W.  In  observations  (Fig.  13a),  the  22-day  band 
propagates  northward  from  40°S  (preeipitation  aetivity 
over  southeastern  subtropieal  South  Ameriea)  to  the 
equator  (Atlantie  ITCZ),  whieh  is  eonsistent  with  the 
results  of  Paegle  et  al.  (2000,  see  their  Fig.  lOf).  In  total, 
6  of  the  14  models  simulate  eoherent  northward  propa¬ 
gation  [i.e.,GFDL2.1,GISS-AOM,MIROC3.2(medres), 
MIROC3.2(hires),  MPI,  and  CSIRO],  but  the  propa¬ 
gation  often  stops  at  10°S,  whieh  is  eonsistent  with  the 
southward  shift  of  the  Atlantie  ITCZ  in  the  models 
(Figs.  1  and  2).  Seven  models  produee  standing  oseillation 
(i.e.,  GFDL2.0,  CCSM3,  PCM.  MRI,  CGCM,  IPSL,  and 
CNRM),  and  one  model  displays  different  propagation 
direetion  in  different  regions  (i.e.,  GISS-ER). 

4.  Summary  and  discussion 

This  study  evaluates  the  intraseasonal  variability  asso- 
eialed  with  the  summer  preeipitation  over  South  Ameriea 


in  14  IPCC  AR4  eoupled  GCMs.  The  results  show  that 
in  the  southern  summer  (November-April),  most  of  the 
models  underestimate  seasonal  mean  preeipitation  over 
eentral-east  Brazil,  northeast  Brazil,  and  the  SACZ. 
Most  models  produee  and  Atlantie  SACZ  to  the  south  of 
that  observed.  Most  of  the  models  eapture  both  the 
40-  and  22-day  bands  around  Uruguay,  but  with  fewer 
aetive  episodes  than  observed.  The  models  also  tend  to 
underestimate  the  total  intraseasonal  (10-90  day)  vari¬ 
anee,  the  40-day  band  varianee,  and  the  22-day  band 
varianee.  In  the  40-day  band,  10  of  the  14  models  sim¬ 
ulate  to  some  extent  the  three-eell  pattern  around  South 
Ameriea,  and  six  models  reproduee  its  teleeonneetion 
with  preeipitation  in  the  south-eentral  Paeifie,  but  only 
one  model  simulates  the  teleeonneetion  with  the  MJO  in 
equatorial  Paeifie,  and  only  three  eapture  its  northward 
propagation  from  50°  to  32°S.  Of  the  seven  models  with 
three-dimensional  data  available,  only  one  reproduees 
well  the  deep  baroelinie  vertieal  strueture  of  the  40-day 
band.  For  the  22-day  band,  only  6  of  the  14  models 
eapture  its  northward  propagation  from  the  SACZ  to 
the  Atlantie  ITCZ. 

Faetors  hypothesized  to  be  important  for  simulating 
subseasonal  variability  inelude  air-sea  interaetion,  land- 
atmosphere  interaetion,  model  resolution,  and  model 
physies.  Regarding  air-sea  interaetion,  all  models  ana¬ 
lyzed  in  this  study  are  eoupled  GCMs,  but  they  still  have 
signifieemt  diffieulties  in  simulating  the  subseasonal  vari¬ 
ability.  However,  previous  studies  have  shown  that  the 
effeets  of  eoupling  depend  strongly  on  the  baekground 
state  (e.g.,  Inness  et  al.  2003;  Turner  et  al.  2005).  Without 
detailed  experimentation  using  eoupled  and  uneoupled 
versions  of  the  same  model  with  similar  mean  states,  few 
firm  eonelusions  ean  be  drawn.  Moreover,  sinee  most 
eoupled  models  are  only  exehanging  air-sea  or  air-land 
fluxes  onee  every  24  h,  more  frequent  eoupling  may  be 
neeessary. 

Land-atmosphere  interaetion  may  also  play  an  im¬ 
portant  role  in  simulating  the  intraseasonal  variability  in 
the  monsoon  regions  (e.g.,  Webster  1983).  In  an  obser¬ 
vational  study,  evidenee  was  found  that  the  40-day  band 
eould  be  loeally  exeited  by  interaetion  with  the  land 
surfaee  states  and  fluxes  in  the  Amazon  rain  forest. 
Future  studies  are  needed  to  assess  how  well  the  IPCC 
models  simulate  the  land-atmosphere  interaetion  over 
the  Amazon  rain  forest. 

Regarding  model  resolution,  we  have  only  one  pair  of 
similar  atmospherie  models  but  with  different  resolution: 
MIROC3.2(hires;  T106)  versus  MIROC3.2(medres;  T42). 
Higher  model  resolution  is  assoeiated  with  weaker  vari¬ 
anee  of  the  40-day  band  (Fig.  6),  but  stronger  varianees  of 
the  22-day  band  (Fig.  12).  It  improves  the  propagation  of 
the  40-day  band  (Fig.  7)  but  not  the  22-day  band  (Fig.  13). 
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Fig  13.  Lag  correlation  of  the  22-day  band  precipitation  anomaly  averaged  between  50“  and  60“W  with  respect  to 
itself  at  30“S,  55“W.  Shading  denotes  the  regions  where  lag  correlation  is  above  the  95%  confidence  level. 
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However,  these  results  may  be  model  dependent,  sinee 
the  resolution  dependenee  is  often  related  to  the  speeifie 
eharaeteristies  of  model  physies.  Moreover,  model  res¬ 
olution  also  affeets  the  representation  of  topography, 
sueh  as  the  Andes  Mountains,  whieh  may  alter  the 
Rossby  wave  trains  that  enter  into  South  Ameriea  from 
the  South  Paeifie.  Sinee  all  of  our  models  are  in  a  rela¬ 
tively  low  resolution  (Table  1),  with  the  highest  resolu¬ 
tion  (T106)  being  about  125  km,  the  poor  representation 
of  the  Andes  Mountains  may  eontribute  to  the  model 
limitations  in  eorreetly  representing  subseasonal  vari¬ 
ability  in  South  Ameriea. 

Regarding  model  physies,  an  interesting  finding 
of  this  study  is  that  the  six  models  simulating  large 
total  intraseasonal,  40-  and  22-day  band  varianees 
[i.e.,  MIROC3.2(hires),  MPI,  CNRM,  GFDL2.0, 
MIROC3.2(medres),  and  CSIRO]  are  just  the  models 
produeing  large  varianees  for  the  eonveetively  eoupled 
equatorial  waves  in  the  tropies  (Lin  et  al.  2006).  A  eom- 
mon  eharaeteristie  of  these  models  is  that  there  is  some 
form  of  moisture  trigger  assoeiated  with  their  eonveetion 
seheme.  We  have  eondueted  a  series  of  GCM  sensitivity 
experiments  to  test  the  effeets  of  moisture  trigger  on 
the  simulated  intraseasonal  variability  assoeiated  with 
summer  preeipitation  over  South  Ameriea  in  the  Seoul 
National  University  GCM.  Three  different  eonveetion 
sehemes  are  used  ineluding  the  simplified  Arakawa- 
Sehubert  (SAS)  seheme,  the  Kuo  (1974)  seheme,  and 
the  moist  eonveetive  adjustment  (MCA)  seheme,  and  a 
moisture  eonveetive  trigger  with  variable  strength  is 
added  to  eaeh  seheme.  The  results  show  that  adding  a 
moisture  trigger  signifieantly  enhanee  the  varianees  of 
both  the  40-  and  22-day  bands.  The  results  will  be  re¬ 
ported  in  a  separate  study. 
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